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We compare topological insulator materials and Rashba coupled surfaces as candidates for engi- 
neering p+ip superconductivity. Specifically, in each type of material we examine 1) the limitations 
to inducing superconductivity by proximity to an ordinary s-wave superconductor, and 2) the ro- 
bustness of the resulting superconductivity against disorder. We find that topological insulators 
have strong advantages in both regards: there are no fundamental barriers to inducing supercon- 
ductivity, and the induced superconductivity is immune to disorder. In contrast, for Rashba coupled 
quantum wires or surface states, the the achievable gap from induced superconductivity is limited 
unless the Rashba coupling is large. Furthermore, for small Rashba coupling the induced super- 
conductivity is strongly susceptible to disorder. These features pose serious difficulties for realizing 
p-|-ip superconductors in semiconductor materials due to their weak spin-orbit coupling, and suggest 
the need to seek alternatives. Some candidate materials are discussed. 



I. INTRODUCTION 

Superconductors with p + ip pairing symmetry have 
long been expected to posses zero-energy Majorana 
bound states in vortex cores^ or at the ends of one- 
dimensional structures^. These Majorana fermion bound 
states are expected to exhibit non-Abelian exchange 
statistics"'^'^, and have been proposed as a basis for topo- 
logical quantum computers which would be protected 
from decoherence"^"^ . Consequently, there is a growing 
interest in realizing a robust p + ip superconductors in 
the laboratory. Such p + ip superconductors are thought 
to naturally occur in triplet paired fermionic superflu- 
ids (such as '^He A or Sr2Ru04)^'^, and in the PfafRan 
quantum Hall state® aX v = 5/2. However, these systems 
are all experimentally delicate, and despite extensive ex- 
perimental work, direct evidence of Majorana fermions 
remains elusive. 

Recently, the possibility of engineering effective p + 
ip superconductors in more conventional materials has 
arisen^^^^. A common thread in these proposals is the 
use of spin-orbit coupling to convert conventional super- 
conductivity into p + ip superconductivity, typically by 
inducing s-wave superconductivity in a 2D material with 
spin-orbit coupling by proximity to an ordinary bulk su- 
perconductor. Among these proposals, two dominant 
classes of candidate materials have emerged: 1) surface 
states of topological insulator (TI) matcrials^^ and 2) 
semiconducting quantum-wires or two-dimensional elec- 
tron gases (2DEGs) with Rashba spin-orbit coupling Ur 
and induced magnetization In this work, we 

provide a detailed comparison of induced superconduc- 
tivity in these two classes of materials and discuss the 
comparative advantages and disadvantages of using each 
of these materials to construct a, p + ip superconductor. 

We first examine the prospects for inducing supercon- 
ductivity in TI surface states and Rashba materials by 
the proximity effect. For TI surface states, the induced 
s-wave pairing is always converted into p + ip pairing due 
to the topologically protected winding of the TI surface 



Bloch wave-functions. Consequently, for a sufficiently 
good interface between the TI surface and a bulk su- 
perconductor, it is possible to induce the full bulk pair- 
ing gap Aq on the TI surface^®. The situation is more 
complicated for Rashba materials, where a more delicate 
balance of spin-orbit and magnetization is required to 
achieve p-l-ip superconductivity. For these materials, the 
size of the induced superconducting gap is limited not 
only by the transparency of the interface to the bulk su- 
perconductor, but also by the magnetization and Rashba 
energy scales which we denote by Vz and Ur respec- 
tively (see Section II B for detailed definitions). In par- 
ticular for small Rashba coupling {Ur <^ T4,Ao), we 
find that the induced superconducting gap is limited to 
iyc7^< Ao. 

We then analyze the effects of disorder on the induced 
superconductivity. Since superconductivity in the 2D 
surface layer is induced by proximity rather than spon- 
taneously developed by phonon interaction, and since 
the induced superconductivity has s-wave symmetry one 
might expect that the disorder cannot reduce the induced 
pairing gap. On the other hand, disorder is pair-breaking 
for p-wave superconductors, and the induced supercon- 
ductivity is effectively converted into p+ip superconduc- 
tivity. Therefore it is not a priori clear what the effect 
of disorder will be. By computing the disorder averaged 
density of states, we find that, for TI materials, the in- 
duced superconductivity is immune to disorder, and ar- 
gue on general grounds that this immunity is a direct 
consequence of time-reversal invariance. In the Rashba 
2DEG's however, the induced magnetization required to 
realize a single helicity p + ip superconductor inherently 
breaks time-reversal symmetry leaving the induced su- 
perconductivity vulnerable to disorder. We find that, 
while the induced superconducting gap A never fully 
closes from disorder, it can be sharply reduced from its 
clean value. The degree of vulnerability to disorder de- 
pends again on the size of the Rashba coupling Ur. For 
small Ur, a is strongly suppressed even for very weak 
disorder for which the superconducting coherence length 
^0 is only a few percent of the mean-free path £. This 
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sharp decrease is more drastic than the case of magnetic 
impurities in a conventional superconductor, for which 
superconductivity is destroyed only when ^ 1. 

In both regards, the TI materials offer advantages, al- 
lowing robust induced superconductivity that is immune 
to disorder. This suggests that TI materials may there- 
fore be the most promising route to realizing topological 
superconductivity and Majorana fermions. However, im- 
plementing a p + ip superconductor using a TI surface 
state requires many further developements in material 
growth and interface engineering, and therefore it may 
still be desirable to work with more conventional mate- 
rials with strong Rashba splitting. 

So far, the theoretical and preliminary experimental 
work on building a, p + ip superconductor from Rashba 
2DEGs has largely focused on semiconductor materials 
and in particular on semiconductor nanowires^^. How- 
ever, in light of our analysis, the very low Rashba en- 
ergy scales in semiconductors raise serious challenges for 
inducing superconductivity. Namely, small Ur greatly 
limits the size of the induced superconducting gap and 
furthermore, renders the resulting superconductor ex- 
tremely sensitive to disorder. These drawbacks suggest 
that an alternative class of materials with stronger spin- 
orbit coupling should be sought. 

Extremely large Rashba splittings on the order of leV 
have been observed in surface alloys of metals and heavy 
elements such as Bi on Ag(lll)^^. In an earlier paper^'^, 
we had proposed this surface alloy as a promising can- 
didate material. However, this extreme Rashba strength 
can also create new problems. Namely, large Rashba cou- 
pling leads to large carrier density, making it difficult to 
adjust the chemical potential by gating. This is problem- 
atic because one must be able to fine-tune the chemical 
potential to achieve topological superconductivity and to 
manipulate Majorana end-states. It is therefore desire- 
able to find materials with strong enough Rashba cou- 
plings to avoid problems with induced superconductivity 
and disorder, but not so strong that gating becomes im- 
possible. 

One particularly promising candidate is the (110) sur- 
face of Au, which first-principles calculations predict will 
exhibit surface bands with sizeable Rashba splitting^°. 
These surface bands naturally lie within < 50meV of the 
bulk Fermi level, indicating that it should be possible to 
move the chemical potential into the topological regime 
using gating. Furthermore, the crystal symmetry of the 
(110) surface allows for a combination of Rashba and 
Dresselhaus type spin-orbit couplings. If both types of 
spin orbit coupling are present, magnetization could be 
induced by applying an in-plane field rather than by prox- 
imity to a ferromagnetic insulator^^. This would greatly 
simplify the proposed setup for realizing Majoranas. 

This paper is organized as follows: we begin with a 
review of the proposed route to engineering a p + ip su- 
perconductor from TI and Rashba 2DEG materials. We 
then introduce a simple model of the proximity effect in 
these materials, and show how to choose system parame- 



ters in order optimize the induced superconducting gap. 
Subsequently, we turn to the issue of disorder, and de- 
rive the disorder averaged Green's functions and density 
of states for weak to moderate disorder {kp£ ^ 1). Fi- 
nally, we close with a discussion of the relative strengths 
and weaknesses of each class of materials for realizing a 
p + ip superconductor. 



II. OVERVIEW OF PROPOSED ROUTES TO 
TOPOLOGICAL SUPERCONDUCTIVITY 

In this section we briefiy review the proposed routes 
to creating a topological superconductor from the sur- 
face state of a bulk topological insulator (TI) or from a 
2DEG with strong Rashba spin-orbit coupling. In both 
classes of materials, spin-orbit coupling creates surface 
bands in which electron spin is locked with respect to 
the direction of propagation. This helical locking of elec- 
tron spin direction to propagation directions causes the 
electron spin to wind as one traverses a loop around the 
Brillouin zone. Consequently, if s-wave superconductiv- 
ity is induced by proximity to a bulk superconductor, the 
helical winding of the surface Bloch wave-functions effec- 
tively converts the induced superonductivity into p + ip 
superconductivity. Specifically, when re-expressed in the 
basis of the surface bands, the induced s-wave pairing 
term takes the form of a p-wave pairing term^^'^^. 

Throughout the paper, we work in the basis of time- 
reversed pairs: 
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which is convenient for discussing superconductivity. 
Here we take the usual representation T = —iaylC of the 
time-reversal operator, where {(7x,y,z} are Pauli matrices 
in the spin basis, and /C denotes complex conjugation. 
We consider inducing the pairing term: 

i?A = A^4_/_,^-f /i.e. 

k 

= J2M'UTi>k)^ + h.c.^Y.'^iAr^^, (2) 



by proximity to an ordinary superconductor, where 
{'''1.2,3} are Pauli matrices in the particle- hole basis. Here 
we have chosen a gauge in which the pairing order pa- 
rameter A is purely real (this is justified, as we are not 
presently concerned with situations where the supercon- 
ducting phase is inhomogenous or fluctuating). 

In proposals to realize Majorana fermions, there are 
two relevant energy scales protecting the coherence of 
information stored among the Majorana fermions. The 
first is the (p-wave component of the) induced supercon- 
ducting pairing gap. A, which sets the energy scale for 
single particle excitations that can change the fermionic 
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parity of a pair of Majoranas. The second is the so- 
caUed mini-gap to locahzed excitations near each Majo- 
rana. Since locahzed excitations cannot change non-local 
fermion parity, the mini-gap is important only when two 
Majoranas are brought close to each other for measure- 
ment purposes^^. For Majoranas realized as end-states, 
the mini-gap also scales linearly with the induced p-wave 
pairing gap^^. Therefore, in order to perform quantum 
coherent manipulations of Majorana fermions, it is im- 
portant to achieve a robust pairing gap, and to work at 
temperatures much lower than this gap. 




A. Topological Insulators 

The low-energy continuum Hamiltonian, H^i 
Efc *I'HTi(A:)*fe of a TI surface is^^: 
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where v is the Dirac cone velocity. The surface eigen- 
states with energies ±v|A:| form the upper and lower 
branches of a single Dirac cone with definite spin-helicity: 
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where 4>k = tan^-'^ {kx/ky). 

When expressed in the c± basis, the pairing term 
(Eq. 2) takes the form of an ideal spinless p + ip 
superconductor^^ 



h.c. 



(5) 



B. Rashba 2DEG 



FIG. 1. (Color online) Band-structure with Rashba coupling, 
magnetization, and induced superconductivity. Three rele- 
vant energy scales are labelled: 2Abg is the energy gap be- 
tween the e+ and e_ bands at the Fermi surface, Aps is the 
induced p-wave pairing gap at the Fermi surface, and 2Vz is 
the induced magnetization gap. 



OM^k) out of the xy-planc modifying the surface eigen- 
states and corresponding dispersions: 
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Re-expressing H/^ in the eigenbasis of both Rashba and 
Zeeman couplings, one finds that, in addition to p ± 
pairing Ap(fc)k^ ~ {ck,±C-k.±) between fermions both 
in band e±, the canting 9m introduces an s-wave pair- 
ing component As(fc) ~ (cfc +c_fc _) between fermions 
c_(- and c_ in bands £_(. and e_ respectively where: 
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The low-energy continuum Hamiltonian, = 
Efe *I^R (^)*fc of a Rashba coupled 2DEG with induced 
magnetization is: 



■HR(k) = [S_k + aRZ ■ (cr X k)] rg + V^cr^ 



(6) 



I 2 , 

where = — is the spin-independent dispersion, 
/i is the chemical potential, Un = 2ma^ is the Rashba 
coupling strength, and Vz is the induced Zeeman splitting 
responsible for the surface magnetization. 

In contrast to the TI case, without breaking time- 
reversal symmetry, there are two helicities present at 
each energy. Consequently, in order to construct a single 
species p + ip superconductor, it is necessary to explicitly 
break time-reversal symmetry, in this case by introducing 
magnetization, Vz, to remove one of these helicities. 

The Rashba coupling an creates two helical bands with 
energies e^^^ = ^feiQ;ij|fc| and spin-wavefunctions c± (as 



and = {ky ±ikx) /k. As discussed in^-^, one has 
a topological superconductor with potential Majorana 
bound states so long as T4 > A, and so long as /z lies 
within the Zeeman gap (|/i| < 14). It is most advanta- 
geous to set fj. = 0, placing the chemical potential in the 
middle of the Zeeman gap (which can be done either by 
electrostatic gating or chemical doping), and we will take 
fi = throughout the remainder of this paper. 

In this system, there are two excitation gap energy- 
scales: the first is the pairing gap at the Fermi surface 
{k = kp) given by: 
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The second is the Zeeman gap at fc = 0, given (for /i = 0) 
for the TI case). Vz cants the helical bands by angle by jVz — A|. The smaller of these two energy scales sets 
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the bulk gap to single-particle excitations which would 
destroy the non-local information stored among Majo- 
rana bound-states. We note that the relative strength 
of the Rashba spin-orbit coupling Ur and the Zeeman 
splitting Vz determines the size of the pairing gap at fc^?. 
For Vz ^ Ur, the pairing gap is only a small fraction of 
originally induced A. If the Zeeman gap closes (i.e. if 
Vj < A), then both helicities are present and the system 
is topologically trivial. 



III. PROXIMITY INDUCED 
SUPERCONDUCTIVITY 

A. Simple Model of Proximity Effect 

In this section, we consider the interface between a 
bulk s-wave superconductor and cither a topological in- 
sulator surface or a Rashba coupled surface state with 
induced magnetization Vz- As a simple model of this in- 
terface, we consider a bulk superconductor described by 
the BCS Hamiltonian: 



Hb = 
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(11) 



coupled to the surface through a clean planar interface 
described by the bulk-surface tunneling term: 
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which conserves momentum ky parallel to the interface, 
and is independent of the transverse momentum k± per- 
pendicular to the interface. Here bj. ^ and cj. ^ are the 
electron creation operators (with momentum k and spin 
a) for the bulk superconductor and surface respectively, 
Eb is the non-superconducting bulk dispersion which we 
will linearize about the chemical potential and Aq is 
the bulk s-wave pairing amplitude. 

Since surface-bulk tunneling conserves in-plane mo- 
mentum, the bulk tunneling density of states (in the 
absence of superconductivity) is given by the one- 
dimensional expression NB{sB{k)) = {deB{k) / dkz)~ . 
Assuming that Nb varies slowly with energy, Hb-s in- 
duces the following self-energy correction to the surface 
Green's function: 
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where 7 ~ Nb{Q)\T\'^ is convenient measure of the 
strength of surface-bulk coupling corresponding to the 
width of the surface resonance that would result from 
^^B-s without bulk-superconductivity (Aq = 0). 

Incorporating Sr into the surface Green's function 
gives: 
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where Zy is the reduced in quasi-particle weight due to 
the bulk-surface hybridization: 



Zr(iw) = 1 
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The quasi-particle weight can be interpreted as the frac- 
tion of time that a propagating electron resides in the 
surface, as opposed to the bulk. The surface-bulk tun- 
neling induces a pairing term Ati in the surface where: 



A = (1 - Zr)Ac 



(16) 



For strong surface-bulk coupling (7 ^ Aq or equivalently 
Zr <C 1) a sizeable fraction of the bulk pairing is induced 
on the surface. 

However, this is not the only effect of the interface. 
From Eq. 14 we see that the surface-bulk coupling renor- 
malizes the surface Hamiltonian, effectively rcscaling the 
coefBcients by a factor of Zy-: 
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The effects of this renormalization are markedly different 
for topological insulators and Rashba 2DEG's, and we 
will consider each case in turn. 



1. Proximity Effect for TI Surface 

Renormalization of the topological insulator surface 
due to coupling to a bulk superconductor simply results 
in rescaling the Fermi velocity vp — >■ vp = Zrvp and 
chemical potential fj. ^ fl = ZrfJ,. However, the nature of 
the induced pairing is independent of vp and |^|, rather, 
it depends only on the helical spin-winding of the sur- 
face Bloch-wavefunctions as one traverses a loop around 
the single Dirac cone in the surface Brillouin zone. Since 
this helical winding is unchanged by the bulk-surface cou- 
pling, the induced pairing symmetry will be preserved 
regardless of 7. 

Therefore, for TI materials there are no funda- 
mental restrictions to pursuing arbitrarily strong cou- 
pling between the TI surface and the nearby bulk- 
superconductor, and consequently it is in principle possi- 
ble to induce the full bulk-gap Aq on the TI surface. This 
advantageous feature of TI surfaces was first pointed out 
in Ref. 18. In the subsequent section, we will see things 
are not so simple for the inducing superconductivity in a 
Rashba coupled 2DEG. 



2. Proximity Effect in a Rashba Coupled 2DEG 

Inducing superconductivity in Rashba coupled surface 
states is a more delicate matter than in TI surfaces. 
Whereas the helical TI surface states are topologically 
guaranteed to convert induced s-wave pairing into effec- 
tive p + ip pairing, constructing effective p + ip pairing 
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FIG. 2. Optimal parameters for producing a large pairing 
gap in a Rashba coupled surface via the proximity effect as 
a function of Rashba coupling strenght Ur. Panels a.-d. re- 
spectively show the optimal excitation gap Ap°^*', Zeeman 
splitting Vz°'^*\ quasi-particle residue ^("p')^ and surface- 
superconductor coupling energy All energies are mea- 
sured with respect to the pairing amplitude Ao of the bulk 
superconductor. 



in a Rashba coupled 2DEG requires a careful balance of 
spin-orbit coupling to create helical winding and magne- 
tization Vz to remove one of the helicities. In particular, 
the Zeeman gap at fc = is given by \Vz — A | , and must 
not close. Consequently, in order to engineer a topo- 
logical superconductor, one must arrange for the surface 
magnetization to exceed the induced pairing amplitude: 
Vz > A- Since V, = ZrK, and A = (1 - Zr)Ao, one 
immediately sees that a balance must be struck to en- 
sure that the bulk-surface interface is sufficiently strong 
to induce pairing, but not so strong that it destroys the 
magnetization gap at k = 0. 

We now turn to a quantitative analysis of the effect of 
proximity induced pairing in a Rashba 2DEG. The goal 
will be to find the optimum set of parameters in order 
to achieve a large p-wave superconducting gap on the 
Rashba coupled surface. There are two excitation gaps, 
the magnetization gap, — A| = \ZrVz~{l — Zr)Ao\, at 
k = and the p-wave pairing gap Aps = {1 — Zy)Afs at 
the Fermi-surface. The smaller of these two energy scales 
sets the minimum excitation energy gap. The former 
decreases with Z-p, whereas the latter increases with Zr, 
indicating that the optimal Z-p is: 

A.. + A„ ^^^^ 



^(opt) _ 



Vz + Afs + Ao 
which corresponds to an optimal excitation gap of: 

. (opt) _ Vz 



^FS 



Vz 



(19) 



■ Afs + Ao 

where Aps is a function of Ur, Vz, and Ao given by Eq. 
9. 



Figure 2a. shows the optimum achievable value of 
Egup as a function of Rashba coupling strength Ur, and 
Figure 2b-d. show the corresponding optimal values of 
Vz, Zr, and 7. In practice, it will likely not be possi- 
ble to fine-tune the interface transparency, 7, between 
the Rashba surface and the adjacent superconducting 
layer, or the induced magnetization Vz- Rather, these 
parameters will be determined by the detailed structure 
of the surface-superconductor and surface-magnetic in- 
sulator interfaces. Therefore, one should view Figure 2a. 
as an upper bound on the practically achievable induced 
pairing gap. Even so, the general trend is clear: for small 
Rashba splitting, Ur, only a small fraction of the bulk 
pairing gap Aq is induced on the surface, whereas for 
large Ur a substantial fraction of Aq is achievable. 

This analysis highlights one of the potentially serious 
drawbacks of using materials with weak spin-orbit cou- 
pling. In the limiting case of small Ur, it is advantageous 
to arrange Zr — 1/2 , large Aq ^ Ur, and Vz = Aq, in 
which case: 



lim E^gap < 
Ur^O 



VUrAo 



(20) 



In particular for semiconductor materials in which typi- 
cal Rashba couplings are typically of the order of 0.2 — 
O.SA'^-^), even after carefully optimizing Vz, 7, and Ag 
only a p-wave pairing gap on the order of 0.1 — OAK 
is achievable. Such small excitation gaps would require 
operating at temperatures much smaller than O.IK in 
order to avoid thermal excitations, which could pose dif- 
ficulties for experiments. Furthermore, as will be shown 
in more detail below, small Ur puts stringent restrictions 
on sample purity, as even small amounts of disorder will 
further suppress induced pairing. In contrast, the situa- 
tion is much more hopefull for materials with strong spin 
orbit couplings. For strong Rashba couplings, one induce 
nearly the full superconducting gap, Ao, given sufficiently 
transparent superconductor-surface interfaces. 



B. Surface Resonances 

So far, we have been implicitly considering an artifi- 
cial interface between a 2D material (either a TI surface 
or Rashba 2DEG) and a different superconducting mate- 
rial. A potentially simpler alternative for realizing p + ip 
superconductivity, is to use the naturally occuring in- 
terface between a bulk-superconductor and its surface. 
This approach would eliminate the need to find compat- 
ible materials to engineer an appropriately transparent 
interface. 

The formalism developed above applies equally well 
in this case. Namely, if electronic states on the surface 
of a bulk metal occur at the same energy and momen- 
tum as bulk states, then the surface states decay into the 
bulk leaving behind broadened resonances. If the bulk 
becomes a superconductor, the surface-bulk coupling in- 
duces superconductivity on the surface. Denoting the 
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width of the surface-resonance (in the absence of bulk- 
superconductivity) by 7, the induced superconductivity 
is again described by Eqs. 14 and 15. It is also pos- 
sible that surface states coexist at the same energy as 
bulk bands, but reside in regions of the Brillouin zone 
for which there are no bulk-states. In this case there is 
no direct tunneling from the surface into the bulk, and 
the surface state would remain sharp state rather than 
broadening into a resonance. Consequently to obtain su- 
perconductivity on the surface, one would need to rely on 
some scattering process (e.g. phonon, electron-electron, 
or disorder scattering) to transfer electrons between sur- 
face and bulk states. 

For natural superconducting metals with strong spin- 
orbit coupling (such as Pb), the electrostatic potential 
created by the material's surface interupts the bulk inver- 
sion symmetry, giving rise to a surface Rashba coupling. 
If the surface Hamiltonian has appropriate combinations 
of Rashba spin-orbit coupling and magnetization (as de- 
scribed above), then the induced surface superconductiv- 
ity will again have effective p + ip pairing symmetry. 

A related approach is possible for topological insulator 
materials, where it has been demonstrated^^ that doping 
can produce superconductivity with transition tempera- 
tures Tc ^ 0.15 — 5.5K. Furthermore, it is common^'^ 
that samples of materials such as Bi2Se3 that are ex- 
pected to be bulk-insulators, are actually metallic. In 
these "topological metals", the topologically protected 
surface states that would appear for a bulk insulator ap- 
pear instead as resonances^'^. In fact a large amount of 
experimental effort is currently focused on finding materi- 
als with genuinely insulating bulks in order to investigate 
surface electron transport. However, for the purpose of 
engineering a.p + ip superconductor, this surface-bulk co- 
existence is actually advantageous, and the combination 
of bulk superconductivity and surface-bulk coupling will 
result in an effective p + ip superconductor at the surface 
of a superconducting topological metal. 

To examine whether p + ip superconductors built from 
surface resonances also exhibit Majorana bound states, 
for example in vortex cores or at the ends of one- 
dimensional magnetic domains, one can write down the 
T-matrix for scattering from a vortex or domain wall 
and look for poles at zero-energy. For a static vor- 
tex or domain wall configuration, the T-matrix at zero- 
energy is constructed from various products of surface 
Green's functions (see Eq.l4) also at zero-energy. Since 
Sr(w ~ 0) = Z{uj ~ 0)Aori, the surface Green's func- 
tion is identical to that of an ideal p + ip superconductor 
with gap A = ZAq. Therefore surface-resonance p + ip 
superconductors will exhibit zero- energy Majorana bound- 
states under exactly the same conditions as the effective 
p + ip superconductor discussed previously. These Ma- 
jorana states are localized to the surface layer and arc 
protected against decaying into bulk states because of 
the bulk superconducting gap. 
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FIG. 3. Panel A shows a diagrammatic representation of Eq. 
23 for the disorder averaged Green's function and self-energy 
respectively. Disorder scattering is represented by dashed line 
originating from an x. For delta- function-correlated impuri- 
ties only multiple scatterings from the same impurity con- 
tribute. Panel B shows an example of a crossed diagram 
(right) that is sub-leading in (kpi)^^ compared to the non- 
crossed diagram with the same number of disorder scatterings 
(left). 



IV. DISORDER 

In this section we show that disorder effects the in- 
duced superconductivity very differently in the TI scheme 
as opposed to the Rashba scheme. To model disorder, we 
consider a random on-site potential 

Hms = Yl ^(^)4,aCr,. (21) 

that has only short-range correlations 

V{r)V{r') = W^d{r - r') (22) 



where W is the disorder strength and (• • •) indicates 
an average over disorder configurations. It is useful 
to parameterize disorder either by the scattering time 
T = l/iV(0)VF^ or the mean free path i = vpT where 
vp is the Fermi velocity of the surface layer and N{0) 
is the surface-density of states. Furthermore, we con- 
sider moderate disorder that is too weak to induce lo- 
calization, specifically that kp£ ^ 1, but make no other 
assumptions on disorder strength. Since non-planar dis- 
order scattering diagrams are sub-leading in (kpi)^^, the 
regime kp£ ^ 1 allows for a controlled expansion for the 
disorder self-energy. 

The disorder averaged self-energy and Green's function 
are related by the following set of coupled equations: 

G{iuj, k) = \_QQ{iuj, ky^ — T,{ioj)\ 

I](iw) = I¥V3^^(iw,fc)T3 (23) 

k 

where the bare (non-disordered Green's function) Qq is 
given by Eq. 14, and incorporates the proximity induced 
superconductivity. 
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A. Time-Reversal Symmetry 

In the subsequent discussion of disorder, the presence 
or absence of time-reversal (TR) symmetry plays a key 
role. We will presently show that when (TR) symme- 
try is present, induced superconductivity is immune to 
the presence of disorder. The proof of this principle is 
most conveniently conducted in the basis of time reversed 
pairs^'' (see Eq. 1), in which the Hamiltonian for the 
disordered system with time-reversal symmetric s-wave 
pairing induced by proximity effect can be written as the 
following block-matrix: 



H = 



Ho + V AI 
AI -{Ho + V) 



(24) 



Here Hq is the Hamiltonian of the (clean) surface, V is 
random on-site disorder, A is the induced pairing am- 
plitude, I is the N X N identity matrix where N is 
the number of degrees of freedom in the system. Since 
AI commutes with Hq and V, the eigenvalue problem 
det {H — e) = can be simplified: 



= det {H -e)^ det (e^I - {Hq + Vf 



(25) 



Denoting the eigenvalues oi Hq + V by {e„}, the eigen- 
values of H are ±£^„ where: 



(26) 



which is bounded below by A, independent of the par- 
ticular disorder configuration. 

These manipulations show that, so long as the 2D sur- 
face Hamiltonian is TR invariant, disorder cannot reduce 
the superconducting gap. As an aside, it is useful to note 
that the above considerations do not depend on V being 
spin-independent so long as it preserves TRI. In partic- 
ular strong spin-orbit impurity scattering will also not 
reduce the superconducting gap. 



B. Disordered Topological Insulators 

Since the TI Hamiltonian (Eq. 3) is TR invariant, from 
the previous discussion we know that the pairing gap can- 
not be diminished by disorder. As a simple demonstra- 
tion of this general principle, one can explicitly calculate 
the disorder averaged Green's function and self-energy 
(see Eq. 23). 

In typical experimental situations, the TI surface 
states are intrinsically doped away from the surface Dirac 
point leaving an appreciable density of states at the Fermi 
surface^'^. In this case, /i ^ A, and one finds: 



S(iu;) 



VA2 



: {iuj — Ati) 



(27) 



where r ^ = 7rA^(0)W^^ is a measure of the disorder 
strength. Incorporating this disorder self-energy into the 



Green's function results in a disorder averaged Green's 
function of the same form as the bare Green's function 
Go, except with a reduced quasi-particle weight Z due to 
disorder scattering: 



g{iuj) = - 



Ziiui) 



iuj — ZUti — Ati 
1 



VA2 



- uj- 



(/i > A) (28) 



For chemical potential tuned to the Dirac point, the re- 
sults are similar except that the quasi-particle weight 
takes a different form: 



Z(iuj) = 



1 



rhi 



A2 



(a' = 0) (29) 



In either case, inspection of Eq. 28 reveals that the min- 
imum excitation gap. A, is unchanged by disorder scat- 
tering, in agreement with the general principles outlined 
above for TR invariant systems. 



C. Disordered Rashba 2DEGs 

In contrast to the TI case, creating a topological su- 
perconductor from a Rashba 2DEG requires explicitly 
breaking TR symmetry by inducing surface magnetiza- 
tion Vz- Without TR symmetry, the general arguments 
outlined above do not apply, and the induced pairing 
is vulnerable to disorder scattering. The analysis be- 
low demonstrates that the pair-breaking effects of disor- 
der scattering are dramatically enhanced by the singular 
density of states at the superconducting gap edge, and 
furthermore that these effects are especially pronounced 
in systems with weak Rashba coupling. Since the pairing 
on the surface is induced by the bulk, it never vanishes. 
However, we shall see that pairing can be greatly reduced 
by even a small amount of disorder unless Ur 3> T4- In 
fact, when Ur <C Vz, the suppression due to disorder is 
more severe than for conventional superconductors with 
magnetic impurities for which the superconducting gap 
typically closes ioi £^o/£ ^ 1. We will see, for Ur ^ Vz, 
that disorder strongly suppresses the induced supercon- 
ductivity even for very weak disorder. 

For superconductivity induced by proximity effect, the 
surface state wave-functions are localized to the surface, 
but extend into the superconductor with characteris- 
tic lengthscale Therefore electrons residing in the 
Rashba material are scattered not only by impurities in 
the Rashba material and interface roughness, but also by 
impurities in the superconductor. For example, even if 
one starts with a pristine semiconductor structure (such 
as a self-assembled nanowire), if superconductivity is in- 
duced by proximity to superconductor with some impu- 
rities then the mean-free path will be set by the super- 
conductor rather than the semiconductor. 

To better understand the distinction between impu- 
rities residing in the Rashba coupled surface and those 
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in the superconductor, we consider each separately. To 
treat either case we find that it is sufficient to replace the 
bare disorder scattering time r^^ = ttN{0)W^ (where 
A''(0) is the surface density of states) by an effective dis- 
order scattering time 



^cff 



(1 - Zrfr-^ 



Surface Disorder , . 
Bulk Disorder ^ ' 



An explicit derivation of these expressions is given in Ap- 
pendix A, but the effective scattering time can be under- 
stood more simply as follows. The fraction of the surface- 
resonance wave-function which lies on the surface is Zr 
whereas the fraction residing in the bulk superconductor 
is (1 — Zr). Therefore, the disorder scattering matrix 
elements should be weighted by either Zr or (1 — Zr) 
for surface and bulk disorder respectively. Eq. 30 can 
then be simply understood by noting that the scattering 
time is proportional to the square of the disorder matrix 
element. By inspection, we see that effects of surface 
disorder are suppressed in the limit of strong surface- 
bulk tunneling, whereas the effects of bulk disorder are 
suppressed in the limit of weak tunneling. Using this ef- 
fective scattering time, we now turn to the problem of 
solving the self-consistency relations given in Eq. 23 us- 
ing the surface Green's function in Eq. 14 which includes 
the effects of proximity to the bulk superconductor. 

For ordinary disordered superconductors, the strength 
of disorder is conveniently parametrized by the ratio of 
the coherence length = ttwf/A to the mean free path 
£ = vpT. For the proximity induced superconductivity 
these parameters are renormalizcd by surface-bulk cou- 
pling and also depend on the type of disorder (surface 
or bulk). We will see that the effects of disorder de- 
pends on disorder strength only through the ratio /^cff , 
where = ttvf/^ is the surface coherence length and 
^cff = vpTcS is the effective mean-free path for disorder 
electrons. This effective ratio can be written in terms 
of the intrinsic ratio of the intrinsic mean-free path, £ 
(un-renormalized by proximity induced superconductiv- 
ity) and coherence length of the bulk superconductor, 
as follows: 



4ff 



^0 



1 - 
(1 



Zr £ 



Surface Disorder 



Bulk Disorder 



By working in terms of this effective ratio, it is possible 
to treat both the cases of surface-impurities and bulk- 
impurtics on equal footing. 



1. Analytic Expressions for Weak Disorder 

For weak disorder {£,o/Ich ^ 1), it is sufficient to eval- 
uate the self-energy to lowest order in disorder scattering 
strength, corresponding to the first diagram for the self- 



energy shown in the top line of Fig. 3: 

eW(*^) = W^'^3 5I^o(»w,fc)r3 



(32) 



Here we emphasize that the value of should be appro- 
priately renormalized according to Eq. 30 depending on 
whether the scattering considered occurs in the surface 
or in the bulk superconductor. For Apg ^ V^, the dom- 
inant contributions to the fc-integral come from near the 
Fermi-surface. Linearizing the Bugoliuobov dispersion 
about the Fermi-surface, and performing the integration 
yields: 



3fieS(i)(iw) 



1 , 14 

2 A2 



^BG 



^BG 



- A2 UrT3 



TiN{Q)W^ 



+ 2[/^Ati] 



'off 



A2 



(33) 
(34) 

(35) 
(36) 



where r^g is a measure of the disorder strength, given 
by Eq. 30. 

This self-energy alters the spectrum of the disorder av- 
eraged BdG Hamiltonian. For weak disorder, we expect 
the gap at the Fermi surface to change only slightly. To 
find the correction to Aps due to disorder, one needs 
to analytically continue the self-energy to real frequency, 
and then look for a pole in the disorder averaged Green's 
function at lu ~ Aps — 6ui, i.e. to solve 



= dct 



6uj-n{kp) (uj 



to leading order in Suj one finds: 

6u} = fjs'i) (uj = Aps - ScSj 



(37) 



where *o = ("t H '"i ~*^Tj 
H^kp) with energy Aj;-^. 

In the limiting case where 3> Ur 



(38) 

is the eigenvector of 
^-0 ~ 



(31) 75 (0 1 1 



Using Aps 




cn 



(39) 



W,\/2Afs6uj 
and solving for 5uj gives the fol 



lowing expression for the disorder renormalized pairing 
gap at the Fermi-surface: 



A 



FS 




1 - 



2/3' 



_R«offy 



(14 » Ur) 

(40) 
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The unusual non-analytic dependence on disorder 
strength stems from the singular behavior of a; as w — >■ 
Aps, which in turn reflects the Van- Hove singularity in 
the superconducting density of states at the gap edge. 
This Van- Hove singularity enhances the effective disorder 
strength x, and in particular leads to an infinite slope of 
^fs{tJ) as r-gi ^0. 

In the opposite limit, where Ur ^ Vz, — 
i (f — f f — l) and consequently the weak disorder 
correction to the gap energy vanishes to leading order. 
Including sub-leading contributions in Ur/Vz results in: 



ttt 



1 - 



2/3 



{Ur » Vz) 
(41) 



2. Numerical Solution for Moderate Disorder 

For stronger disorder, Eq. 23 must be solved self- 
consistently, which can be done numerically. In order 
to regulate the numerical integrals in the UV we replace 
the continuum dispersion with a periodic one of the form 
^fc = —2t cos(fc) which naturally introduces a finite band- 
width. The top and bottom panels of Figure 4 show the 
dependence of the induced superconducting gap on dis- 
order strength for Ur > Vz and Vz ^ Ur respectively. 

For very weak disorder, ^ ^off: the excitation gap 
exhibits non-analytic infinite initial slope predicted by 
Equations 40 and 41. Stronger disorder never fully closes 
the superconducting gap, however for Vz 3> Ur, the gap 
is largely suppressed even when is only a few percent 
of £cS- In most cases, i?gap is suppressed smoothly with 
increasing disorder strength, however, the £'gap/A curves 
for Ur ~ Vz have a knee-shaped kink at ^o/^cff — 0.07 
after which -Egap drops abruptly. This knee occurs when 
disorder reduces of the magnetization gap at fc = below 
the pairing gap Aps at the Fermi-surface. 

These results can be readily applied both to the case 
where disorder scattering is due to the adjacent super- 
conductor, and when disorder scattering occurs in the 
surface material, by choosing the form for ^o/£eS from 
Eq. 31. It is possible to reduce the sensitivity to bulk- 
disorder by reducing the surface-bulk tunneling strength, 
r. However, reducing F will also cause a smaller prox- 
imity induced pairing gap A. Similarly, it is possible to 
reduce the sensitivity to surface disorder by increasing 
the surface-bulk tunneling rate, though, doing so will be 
detrimental if the surface layer is cleaner than the bulk 
superconductor . 

Here we see a second drawback of using materials with 
low Rashba coupling: in addition to limiting the size of 
the induced pairing gap in the absence of disorder, small 
Rashba coupling renders the topological superconductor 
susceptible even to small amounts of disorder (^o /^cff ^ 
1). While bulk semiconductors arc typically cleaner than 



u rv 

R 



V rv =1/4 

R z 



0.08 




FIG. 4. The excitation gap -Bgap as a function of coher- 
ence length ^0 ~ TTVp/A to the effective mean free path 
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vpT^ff ■ Egap is obtained from numerically solving 



Eq. 23 for a Rashba 2DEG with induced magnetization 
Vz and superconductivity A. Here (• ■ •) indicates renormal- 
ization due to the proximity effect. The effective disorder 
strength has a different form depending on whether disorder 
scattering occurs predominantly in the surface-layer or in the 
bulk superconductor. Both cases can be treated by choos- 
ing the appropriate expression for S^o/lcS from Eq. 31. The 
parameters used in this simulation were t = 1, Vz = 0.1, 
A = 0.01, and various values of Ur. The top panel shows 
curves for Vz S> Ur, the regime appropriate for semiconduc- 
tor materials, whereas the bottom panel shows curves in the 
Ur > Vz regime which could be achieved by using metallic 
thin films with stronger spin-orbit coupling. The magnetiza- 
tion Vz breaks time reversal symmetry rendering the induced 
pairing susceptible to disorder. For 14 S> Ur the gap is al- 
ready strongly supressed when is only a few percent of IcS- 



metallic thin films, their extreme sensitivity to disorder 
will likely be problematic. In particular, great care would 
need to be taken to limit interfacial roughness between 
the semiconductor and adjacent bulk superconductor and 
magnetic insulating film. 

Before concluding, we remark on two possible exten- 
sions of this analysis. Firstly, the effects of disorder were 
treated for fully two-dimensional structures, whereas Ma- 
jorana fcrmions emerge in one-dimcnsional (or quasi-one- 
dimensional) geometries. The effects of disorder in quasi- 
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one-dimensional Rashba coupled structures were ana- 
lyzed numerically in Ref. 13 and 14, and give similar 
results to those given above for two-dimensions. Finally, 
while this analysis has been carried out for the case of 
Rashba-type spin-orbit coupling, we expect similar re- 
sults for systems in which both Rashba and Dresselhaus- 
typc spin-orbit couplings are present. The relevant factor 
in either case is the presence of magnetization Vz which 
breaks time-reversal symmetry and renders the induced 
superconductivity susceptible to disorder regardless of 
spin-orbit type. 



V. DISCUSSION AND CONCLUSION 

In conclusion, we have compared the prospects for con- 
structing an effective p + ip superconductor from TI and 
Rashba 2DEG based materials. We have focused on tech- 
nical limitations to inducing superconductivity in these 
materials, and examined the effects of disorder on the 
induced superconductivity. In both regards, the TI ma- 
terials offer natural advantages. In particular, the effec- 
tive p + ip nature of induced superconductivity in a TI 
surface is guaranteed by the intrinsic spin-helicity of the 
bare TI surface states. As a consequence, there are no 
fundamental limitations for inducing superconductivity 
by the proximity effect. Furthermore, since time-reversal 
symmetry remains intact in the TI surface, the induced 
superconductivity is guaranteed, based on general prin- 
ciples, to be immune to disorder. 

While TI surface states offer certain advantages, TI 
materials are relatively new and many materials chal- 
lenges remain. It may therefore be desireable to construct 
ap+ip superconductor from more conventional materials 
with strong Rashba spin-orbit coupling. Here one needs 
to induce superconductivity by proximity to a conven- 
tional superconductor, and to induce magnetization for 
example by proximity to a ferromagnetic insulator. In 
this case one must strike a comparatively delicate bal- 
ance of spin-orbit coupling and induced magnetization 
to ensure that the resulting superconductor is effectively 
p + ip^^. 

In this regard, materials with small Rashba coupling 
strengths face serious difficulties: 1) the size of the in- 
duced superconducting gap is limited by the size of the 
Rashba coupling and 2) for weak Rashba coupling the 
induced superconductivity is quite fragile and is strongly 
suppressed even by small amounts of disorder. Our anal- 
ysis indicates that an alternative class of materials with 
stronger spin-orbit coupling should be sought. For exam- 
ple, metallic thin films with heavy atomic elements can 
give orders of magnitude larger spin-orbit couplings than 
the semiconductor materials that have so far dominated 
the theoretical discussion. 

Constructing a p + ip superconductor from a Rashba 
2DEG in this way requires engineering a complicated set 
of material interfaces between the Rashba 2DEG, super- 
conductor, ferromagnetic insulator, and gate electrodes. 



However, a fortuitous choice of material could obviate 
the need for the superconducting and ferromagnetic in- 
terfaces. For example, one could try to find bulk metals 
with strong bulk spin-orbit coupling that naturally super- 
conduct and possess surface resonances. In such a sce- 
nario, the coupling between surface-and bulk will then 
automatically induce superconductivity on the surface, 
eliminating the need to build an artificial interface for this 
purpose. A further simplification is possible the metallic 
surface which has the appropriate symmetry such that 
both Rashba and Dresselhaus type spin-orbit coupling is 
present^^. In this case, one could induce magnetization 
by applying an external field rather than by depositing 
a ferromagnetic insulator^^. If a naturally superconduct- 
ing material can be found with the appropriate spin-orbit 
coupled surface resonances, this approach might offer the 
simplest route to an artificial p + ip superconductor and 
Majorana fermions. 
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Appendix A: Surface vs. Bulk Disorder 

In the main text, we argued on conceptual grounds, 
that the above analysis for disorder is easily modified 
to separately treat the two distinct cases where disorder 
scattering occurs predominantly in the surface layer or 
the adjacent bulk superconductor, by replacing the bare 
disorder scattering time by an effective scattering time 
weighted by Zp or (1 — Zp)^ respectively. Here we con- 
sider each case separately and provide an explicit demon- 
stration of this claim. 



1. Surface Disorder 

Consider first the case of a pristine superconductor 
so that impurity scattering occurs only in the surface 
layer. In this case the disorder scattering matrix ele- 
ments are proportional to the fraction of the electron 
wave-functions that resides on the surface. More pre- 
cisely, we can calculate self-energy from disorder scatter- 
ing in Eq. 23 using the surface Green's function in Eq. 
14. Equivalently, in Fig. 3, one should take single solid 
lines to be the surface Green's function in Eq. 14. By 
inspection, we see that the effective disorder matrix ele- 
ments are reduced by a factor oi Zr < I. Therefore, the 
disorder self-energy for surface-disorder is suppressed by 
a factor of Zp. 



11 



The analysis in the main text can be easily modified 
to treat the case where disorder occurs predominantly in 
the surface layer by replacing ri for the surface (without 
surface-bulk tunneling) by the effective ratio: 



(Al) 



and keeping the same ratio Ur/Vz (since this ratio is 
unaffected by the surface-bulk coupling). 



2. Bulk Disorder 

When the dominant source of disorder is in the ad- 
jacent bulk-superconductor, the surface electrons must 
first tunnel into the bulk in order to scatter from the 
disorder potential. This leads to an rcnormalized effec- 
tive disorder strength which is different than the bulk 
value. To demonstrate this, we focus on the limit where 
surface-bulk tunneling is much stronger than bulk disor- 
der (7 <^ W), in which case there is typically no more 
than one disorder scattering event per surface-bulk tun- 
neling event. 

In this limit, it suffices to compute Eq. 23 using the 
surface-Green's function in Eq. 14, with an effective dis- 
order vertex is given by the diagram shown in Fig. 5. In 
this figure the circled F indicates surface-bulk tunneling 
and the dashed line indicates scattering from disorder. 
Written in terms of the bulk Green's functions the effec- 
tive disorder vertex for surface states is: 



F(zcj,k||,Q)r3 



■ ^(fe,,k||)T3 - Ati iuj - £,u+qT3 - Ati 

(A2) 

We are interested primarily in low-frequency behavior for 
which the external momenta {kj|,k|| + Q||} lie near the 
surface Fermi- level. For a given Q|| connecting two points 
on the surface Fermi-surface, there are two values of Qz 
for which k and k -I- Q also lie on the surface. Since the 
dominant contributions are for momenta lying near the 



Fermi-surface, we linearize the bulk-dispersion in the z- 
direction, and denote the tunneling density of states by 

/a£B(k)V ^ 



Nb{0) 



. With these approximations, the 



dkz J 

effective disorder potential for surface electrons is inde 
pendent of momentum transfer Qn: 



V{iio) 



7rA^B(0)|r|- 

VA2 + LO^ 



(A3) 



The disorder self-energy given by the diagram in Fig. 
3a. can then be computed using the surface Green's func- 
tion from Eq.23 and the effective disorder strength in Eq. 
A3. When this self-energy is incorporated into the sur- 
face Green's function, it comes with a factor of {Zr)^. 
Since we are interested in the behavior for co <^ A, we 
may neglect the frequency dependence of various quan- 
tities, and we find that the effective disorder scattering 
time for surface-electrons is weighted (compared to the 
bulk quantity) by a factor: 



/cff 



^A^B(o)|r|' 

A + nNB iO)\ 
{l-Zr?(r-^- 



'bulk 



'bulk 



(A4) 



For a transparent surface-bulk interface (7rA^B(0)|rp ^ 
A), this factor approaches unity, justifying the claim 
made in the main text. 



X 

Q 



(kz ,k|^ k+Q 



FIG. 5. Effective disorder vertex for surface-electrons in 
the case where disorder scattering occurs predominantly in 
the adjacent bulk-superconductor. The solid lines are bulk 
Green's functions, circled Fs indicate surface-bulk tunneling 
events, and the dashed hne ending in x indicates disorder 
scattering that transfers momentum Q 
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